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ABSTRACT 

We have developed a technique to systematically discover and study high-redshift 
supernovae that can be used to measure the cosmological parameters. We report 
here results based on the initial seven of >28 supernovae discovered to date in 
the high-redshift supernova search of the Supernova Cosmology Project. We find 
an observational dispersion in peak magnitudes of ctm b = 0.27; this dispersion 
narrows to 0"m s , con = 0.19 after "correcting" the magnitudes using the light-curve 
"width-luminosity" relation found for nearby (z < 0.1) type la supernovae from the 
Calan/Tololo survey (Hamuy et al. 1996). Comparing lightcurve- width-corrected 
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magnitudes as a function of redshift of our distant (z = 0.35-0.46) supernovae 
to those of nearby type la supernovae yields a global measurement of the mass 
density, f^M = 0.88 IIq'qq for a A = cosmology. For a spatially flat universe (i.e., 

+ = 1), we find Qyi = 0.94 1^28 or ' equivalently, a measurement of the 
cosmological constant, £l\ = 0.06 ^[il (<0.51 at the 95% confidence level). For 
the more general Friedmann-Lemaitre cosmologies with independent Qm and 
the results are presented as a confidence region on the Hm _ ^a plane. This region 
does not correspond to a unique value of the deceleration parameter qo. We present 
analyses and checks for statistical and systematic errors, and also show that our results 
do not depend on the specifics of the width-luminosity correction. The results for 
f^A-versus-fiM are inconsistent with A-dominated, low density, flat cosmologies that 
have been proposed to reconcile the ages of globular cluster stars with higher Hubble 
constant values. 

Subject headings: cosmology: observations — distance scale — supernovae:general 

1. Introduction 

The classical magnitude-redshift diagram for a distant standard candle remains perhaps 
the most direct approach for measuring the cosmological parameters that determine the fate 
of the cosmic expansion (Sandage 1961, 1989). The first standard candles used in such studies 
were first-ranked cluster galaxies (Gunn & Oke 1975, Kristian, Sandage, & Westphal 1978) 
and the characteristic magnitude of the cluster galaxy luminosity function (Abell 1972). More 
recent measurements have used powerful radio galaxies at higher redshifts (Lilly & Longair 1984, 
Rawlings et al. 1993). Both the early programs (reviewed by Tammann 1983) and the recent 
work have proven particularly important for the understanding of galactic evolution, but are 
correspondingly more difficult to interpret as measurements of cosmological parameters. The type 
la supernovae (SNe la), the brightest, most homogeneous class of supernovae, offer an attractive 
alternative candle, and have features that address this evolution problem. Each supernova 
explosion emits a rich stream of information describing the event, which we observe in the form 
of multi-color light curves and time-varying spectra. Supernovae at high redshifts, unlike galaxies, 
are events rather than objects, and their detailed temporal behavior can thus be studied on an 
individual basis for signs of evolution relative to nearby examples. 

The disadvantages of using supernovae are also obvious: they are rare, transient events that 
occur at unpredictable times, and are therefore unlikely candidates for the scheduled observations 
necessary on the largest telescopes. The single previously identified high-redshift [z = 0.31) SN la, 
discovered by a 2-year Danish/ESO search in Chile, was found (at an unpredictable time) several 
weeks after it had already passed its peak luminosity (N0rgaard-Nielsen et al. 1989). 
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To make high-redshift supernovae a more practical "cosmological tool," the Supernova 
Cosmology Project has developed a technique over the past several years that allows the discovery 
of high-redshift SNe la in groups of ten or more at one time (Perlmutter et al. 1997a). These 
"batch" discoveries are scheduled for a particular night, or nights, thus also allowing follow-up 
spectroscopy and photometry on the large-aperture telescopes to be scheduled. Moreover, the 
supernova discoveries are generally selected to be on the rising part of the light curves, and can be 
chosen to occur just before new moon for optimal observing conditions at maximum light. 

Since our demonstration of this technique with the discovery of SN 1992bi at z = 0.458 
(Perlmutter et al. 1994, 1995a), we have now discovered more than 28 SNe, most in two batches 
of ~10 (Perlmutter et al. 1995b, 1997b). Almost all are SNe la detected before maximum light 
in the redshift range z = 0.35-0.65. We have followed all of these supernovae with photometry 
and almost all with spectroscopy, usually at the Keck 10-meter telescope. Other groups have 
now begun high-redshift searches; in particular, the search of Schmidt et al. (1997) has recently 
reported the discovery of high-redshift supernovae (Garnavich et al. 1996a,b). 

We report here the measurements of the cosmological parameters from the initial seven 
supernovae discovered at redshifts z > 0.35. Since this is a first measurement using this technique, 
we present some detail to define terms, demonstrate cross-checks of the measurement, and outline 
the directions for future refinements. Section 2 reviews the basic equations of the technique and 
defines useful variables that are independent of Hq. Section 3 discusses the current understanding 
of type la supernovae as calibrated standard candles, based on low-redshift supernova studies. 
Section 4 describes the high-redshift supernova data set. Section 5 presents several different 
analysis approaches all yielding essentially the same results. Section 6 lists checks for systematic 
errors, and shows that none of these sources of error will significantly change the current results. 

In conclusion (Section 7), we find that the alternative analyses and cross-checks for systematic 
error all provide confidence in this relatively simple measurement, a magnitude versus redshift, 
that gives an independent measurement of Qyi and comparable to or better than previous 
measurements and limits. Other current and forthcoming papers discuss further scientific results 
from this data set and provide catalogs of light curves and spectra: Pain et al. (1996) present first 
evidence that high-redshift SN la rates are comparable to low-redshift rates, Kim et al. (1996) 
discuss implications for the Hubble constant, and Goldhaber et al. (1996) present evidence for 
time dilation of events at high redshift. 

2. Measurement of S7a versus Qm from m—z relation 

The classical magnitude-redshift test takes advantage of the sensitivity of the apparent 
magnitude-redshift relation to the cosmological model. Within Friedmann-Lemaitre cosmological 
models, the apparent bolometric magnitude m(z) of a standard candle (absolute bolometric 
magnitude M) at a given redshift is a function of both the cosmological-constant energy density 



- 4 - 



fi A = A/(3H§) and the mass density Om : 

m(z) = M + 5\ogd L (z;n M ,^A,H ) + 25 

= M + 51ogP L (z;0 M ,^A)-51og^o + 25, (1) 

where di is the luminosity distance and T>l = H^dL is the part of the luminosity distance 
expression 14 that remains after multiplying out the dependence on the Hubble constant (expressed 
here in units of km s _1 Mpc -1 ). In the low redshift limit, Equation |l] reduces to the usual linear 
Hubble relation between m and log cz: 

m(z) = M + 5 log cz — 5 log Hq + 25 

= M + 5 log cz , (2) 

where we have expressed the intercept of the Hubble line as the magnitude "zero point" 
A4 = M — 51og-ffo + 25. This quantity can be measured from the apparent magnitude and 
redshift of low-redshift examples of the standard candle, without knowing Hq. Note that the 
dispersions of M. and M are the same, o_m = &m, since 51og-£/o is constant. (In this paper, we use 
script letters to represent variables that are independent of Hq; the measurement of Hq requires 
extra information, the absolute distance to one of the standard candles, that we do not need for 
the measurement of the other cosmological parameters.) 

Thus, with a set of apparent magnitude and redshift measurements (m(z)) for high-redshift 
candles, and a similar set of low-redshift measurements to determine A4, we can find the best fit 
values of S7 A and SIm to solve the equation 

m(z) -M = 5logV L (z;Q M ,Q A ) . (4) 

(An equivalent procedure would be to fit the low- and high-redshift measurements simultaneously 
to Equation |3|, leaving A4 free as a fitting parameter.) For candles at a given redshift, this fit 
yields a confidence region that appears as a diagonal strip on the plane of versus £Im- Goobar 
& Perlmutter (1995) emphasized that the integrand for the luminosity distance T>l depends on 
J7m and f^A with different functions of redshift, 13 so that the slope of the confidence-region strip 
increases with redshift. This change in slope with redshift makes possible a future measurement 



We reproduce the equation for the luminosity distance, both for completeness and to correct a typographical 
error in Goobar & Perlmutter (1995): 

d L (z;n M ,n A ,H ) = c(1 +fl s( VVl / [{l + z') 2 {l + n M z')-z'(2 + z')n A ]~^ dz') , (3) 
Hoa/\k\ V Jo / 

where, for Qm + > 1, S(x) is defined as sin(:r) and k = 1 — SIm — ^a; for Qm + < 1, S(x) — sinh(a;) and k 
as above; and for Qm + Qa = 1, S(x) = x and k = 1. The greater-than and less-than signs were interchanged in the 
definition of S(x) in the printed version of Goobar & Perlmutter, although all calculations were performed with the 
correct expression. 
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of both and separately, using supernovae at a range of redshifts from z = 0.5 to 1.0; see 
Figure 1 of Goobar & Perlmutter. 

Traditionally, the magnitude-redshift relation for a standard candle has been interpreted 
as a measurement of the deceleration parameter, qo, primarily in the special case of a A = 
cosmology where go and J7m are equivalent parameterizations of the model. However, in the 
most general case, qo is a poor description of the measurement, since T>l is a function of 
and Oa independently, not simply the combination qo = Om/2 — Thus the slope of the 
confidence-region strip is not parallel to contours of constant qo, except at redshifts z « 1. We 
therefore recommend that for cosmologies with a non-zero cosmological constant, qo not be used 
by itself to describe the measurements of the cosmological parameters from the magnitude-redshift 
relation at high redshifts, since it will lead to confusion in the literature. 

Steinhardt (1996) recently pointed out that cosmological models can be constructed with 
other forms of energy density besides and such as the energy density due to topological 
defects. These energy density terms will not in general lead to the same functional dependence of 
luminosity distance on redshift. In this current paper, we do not address these cosmologies with 
additional (or different) energy density terms, since this first set of high-redshift supernovae span 
a relatively narrow range of redshifts. Although constraints on these cosmological models can be 
found from this limited data, our upcoming larger data sets with a larger redshift range will be 
much more appropriate for this purpose. 



3. Low-Redshift SNe la and Calibrated Magnitudes 

To measure the magnitude "zero point" A4, it is important to use low-redshift supernovae that 
are far enough into the Hubble flow that their peculiar velocities are not an appreciable contributor 
to the redshift. It is also better if the sample of low-redshift supernovae were discovered in a 
systematic search, since this more closely approximates the high-redshift sample; our high-redshift 
search technique yields a more uniform (and measurable) magnitude limit than typical of most 
serendipitous supernova discoveries. 

The Calan/Tololo supernova search has discovered and followed a sample of 29 supernovae in 
the redshift range z = 0.01-0.10 (Hamuy et al. 1995, 1996). Of these, 18 were discovered within 
5 days of maximum light or sooner. This subsample is the best to use for determining M, since 
there is little or no extrapolation in the measurements of the peak apparent magnitude or the light 
curve decline rate. The absolute .B-magnitude distribution of these 18 Calan/Tololo supernovae 
exhibits a relatively narrow RMS dispersion, o^™ uy = 0.26 mag, with a mean magnitude zero 
point of M B = -3.17 ±0.03 (Hamuy et al. 1996). 

Recent work on samples of SNe la at redshifts z < 0.1 has focussed attention on examples 
of differences within the SN la class, with observed deviations in luminosity at maximum light, 
color, light curve width, and spectrum (e.g. Filippenko et al. 1992a,b; Phillips et al. 1992; 
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Leibundgut et al. 1993; Hamuy et al. 1994). It appears that these deviations are generally highly 
correlated, and — perhaps surprisingly — approximately define a single-parameter supernova family, 
presumably of different explosion strengths, that may be characterized by any of these correlated 
observables. It thus seems possible to predict, and hence correct, a deviation in luminosity using 
such indicators as light curve width (Phillips 1993; Hamuy et al. 1995; Riess, Press, & Kirshner 

1995) , U—B color (Branch, Nugent, & Fisher 1997), or spectral feature ratios (Nugent et al. 

1996) . 



3.1. Light Curve Width Calibration 

The dispersion o"^^° uy = 0.26 mag thus can be improved by "calibrating," using the 
correlation between the time scale of the supernova light curve and the peak luminosity of the 
supernova. The correlation is in the sense that broader, slower light curves are brighter while 
narrower, faster light curves are fainter. Phillips (1993) proposed a simple linear relationship 
between the decline rate Amis, the magnitude change in the first 15 days past B maximum, 
and the peak B absolute magnitude. (In practice, low-redshift supernova light curves are not 
always observed during these 15 days, and therefore the photometry data are fit to a series of 
alternative template SN la light curves that span a range of decline rates; Amis is actually 
found by interpolating between the Amis values of the templates that fit best.) Hamuy et al. 
(1995, 1996) have now fitted a linear relation for all of the 18 Calan/Tololo supernovae that were 
discovered near maximum brightness, and for the observed range of Amis between 0.8 and 1.75 
mag they obtain: 



M B ,covv = (0.86 ± 0.21)(Ami 5 - 1.1) - (3.32 ± 0.05) . (5) 

This fit provides a prescription for "correcting" magnitudes to make them comparable to 
an arbitrary "standard" SN la light curve of width Amis = 1-1 mag: Add the correction term 
AcorV = (—0.86 ± 0.21)(Ami5 — 1.1) to the measured B magnitude, so that m_B )COXr = ms + A^V • 
(We use the {1.1} superscript as a reminder that this correction term is defined for the arbitrary 
choice of light curve width, Amis = 1-1 mag.) The residual magnitude dispersion after adding this 
correction to the Calan/Tololo supernova magnitudes drops from o"^" uy = 0.26 to o^™orr = 
magnitudes. It is important to notice that the magnitude zero point, M.b = —3.17 ± 0.03, 
calculated from the uncorrected magnitudes is not the same as 

•^B 'corr = ~ 3 - 32 ± °- 05 > the 
intercept of Equation || at Amis = 1.1 mag; this simply reflects the fact that Amis = 1-1 is not 

the value for the average SN la. 

Riess, Press, &; Kirshner (1995) have presented a different analysis of this light curve 
width-luminosity correlation that adds or subtracts different amounts of a "correction template" 
to a standard light-curve template (Leibundgut et al. 1991), creating a similar family of broader 
and narrower light curves. They use a simple linear relationship between the amount of this 
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correction template added and the absolute magnitude of the supernova, resulting in a similarly 
small dispersion in the B and V absolute magnitude after correction, o^bv ~ 0.20 mag. More 
recent results of Riess, Press, & Kirshner (1996) show even smaller dispersion (o"]^ 1 ^ ~ 0.12 mag) 
if multiple color light curves are used and extinction terms included in the fit. 

There remains some question concerning the details of the light curve width-luminosity 
relationship. It is not clear that a straight line is the "true" model relating Amis to Mb, nor that 
a linear addition or subtraction of a Riess et al. correction template best characterizes the range of 
light curves in all bands. However, a simple inspection of the absolute magnitude as a function of 
Amis from Hamuy et al. (1995, 1996) shows primarily a narrow Mb dispersion (ctm b ~ 0.2 mag) 
for most of the supernovae, those with light curve widths near that of the Leibundgut standard 
template (Amis = 1-1 ma g)> as we h as a few slightly brighter, broader supernovae and a tail of 
fainter, narrower supernovae. For the purposes of this paper, a simple linear fit appears to be 
sufficient, since the differences from a more elaborate fit are well within the photometry errors. 

To make our results robust with respect to this correction, we have analyzed the data (a) 
as measured, i.e., with no correction for the width-luminosity relation, adopting the uncorrected 
"zero point" of Hamuy et al, Mb = —3.17 ± 0.03; and (b) with the correction and zero point of 
Equation [5] for the five supernovae that can be corrected with the Hamuy et al. calibration. We 
have also compared the "correction template" approach for the supernova for which photometry 
was obtained in bands suitable for this method, following the prescription and template light 
curves of Riess et al. (1996). 



3.2. Stretch-Factor Parameterization 

For the analysis of Amis m this paper, we use a third parameterization of the light curve 
width/shape, a stretch factor s that linearly broadens or narrows the rest-frame timescale of an 
average (e.g. Leibundgut et al. 1991) template light curve. This stretch factor was proposed 
(Perlmutter 1997a) as a simple heuristic alternative to using a family of light curve templates, 
since it fits almost all supernova light curves with a dispersion of < 0.05 mag at any given time in 
the light curve during the best-measured period from 10 days before to 80 days after maximum 
light. (Physically, the stretch factor may reflect a temperature-dependent variation in opacities 
and hence the diffusion time-scale of the supernova atmosphere; see Khohklov, Miiller, & Hoflich 
1993.) 

The stretch factor s can be translated to a corresponding Amis via the best-fit line 

Ami 5 = (1.96 ± 0.17)(s _1 - 1) + 1.07 . (6) 

Using this relation, the best-fit s-factor for the template supernovae used by Hamuy et al. 
reproduces their Amis values within ±0.01 mag for the range 0.8 < Amis < 1-75 mag covered by 
the 18 low-redshift supernovae. This provides a simple route to interpolating Amis f° r supernovae 
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that fall between these templates. 

In our analysis, we use Equation || together with the Hamuy et al. width-luminosity relation 
(Equation |J) to calculate the magnitude correction term, A|^, from the stretch factor, s. Note 
that Equation || is based on the Amis interpolations of Hamuy et al, not on a direct measurement 
of s, and it will be recalculated once the Hamuy et al. light curves are published. However, 
for available light curves we get close agreement (within approximately 0.04 mag) between 
published Amis values interpolated between light curve templates by Hamuy et al. and the values 
interpolated using s and Equation ||. The uncertainty introduced by this translation is much 
smaller than the uncertainties in the measurement of s for the high-redshift supernovae and the 
uncertainties from approximating the width-luminosity relation as a straight line in Equation [|. 

3.3. Color and Spectroscopic Feature Calibrators 

In addition to these parameterizations of the light-curve width or shape, several other 
observable features appear to be correlated with the absolute magnitude of the supernova. 
Vaughan et al. (1995) suggested that a color restriction B — V < 0.25 mag eliminated the 
subluminous supernovae from a sample of low-redshift supernovae, and Vaughan, Branch, & 
Perlmutter (1996) confirmed this with a more recent sample of supernovae. Branch et al. (1996) 
presented a potentially stronger correlation with U — B color, and showed a very small dispersion 
in graphs of Mb or Amis versus U — B. This result is consistent with the variation in UV flux for 
a series of spectra at maximum light, ranging from the broad, bright SN 1991T to the fast, faint 
SN 1991bg, presented in Figure 1 of Nugent et al. (1996). This figure also showed a correlation of 
absolute magnitude with ratios of spectral features on either side of the Ca II H and K absorption 
trough at 3800 A and with ratios of Si II absorption features at 5800 A and 6150 A. 

Such multiple correlations with absolute magnitude can provide alternative methods for 
calibrating the SN la candle. In the current analysis, we use them as cross-checks for the 
width-luminosity calibration when the data set is sufficiently complete. For future data sets they 
may provide better, or more accessible, primary methods of magnitude calibration. 

3.4. Correlations with Host Galaxy Properties 

There have been some indications that the low luminosity members of the single-parameter 
SN la family are preferentially found in spheroidal galaxies (Hamuy et al. 1995) or that the more 
luminous SNe la prefer late spirals (Branch, Romanishin, & Baron 1996). If correct, this suggests 
that calibration within the SN la family, whether by light curve shape, color, or spectral features, 
is particularly important when comparing SNe la from a potentially evolving mix of host galaxy 
types. 
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4. The High-Redshift Supernova Data Set 

4.1. Discovery and Classification 

The seven supernovae discussed in this paper were discovered during 1992-94 in coordinated 
search programs at the Isaac Newton 2.5 m telescope (INT) on La Palma and the Kitt Peak 
4 m telescope, with follow-up photometry and spectroscopy at multiple telescopes, including 
the William Herschel 4 m, the Kitt Peak 2.3 m, the Nordic Optical 2.5 m, the Siding Springs 
2.3 m, and the Keck 10 m telescopes. The light-curve data were primarily obtained in the 
Johnson-Cousins ii-band (Harris set; see Massey et al. 1996), with some additional data points in 
the Mould /, Mould R, and Harris B bands. SN 1994G was observed over the peak of the light 
curve in the Mould / band (Harris set). All of the supernovae were followed for more than a year 
past maximum brightness so that the host galaxy light within the supernova seeing disk could be 
measured and subtracted from the supernova photometry measurements. Spectra were obtained 
for each of the supernova's host galaxies, and for the supernova itself in the case of SN 1994F, 
SN 1994G, and SN 1994an. The redshifts were measured from host-galaxy spectral features, and 
their uncertainties are all < 0.001. Table g lists the primary observational data obtained from the 
light-curve photometry and spectroscopy. 

We consider several types of observational evidence in classifying a supernova as a type la: (a) 
Ideally, a spectrum of the supernova is available and it matches the spectrum of a low-redshift type 
la supernova observed the appropriate number of rest-frame days past its light-curve maximum 
(e.g., Filippenko 1991.) This will usually differentiate types lb, Ic, and II from type la. For 
example, near maximum light SNe la usually develop strong, broad features, while the spectra of 
SNe II are more featureless. Distinctive features such as a trough at 6150A ( now believed to be 
due to blueshifted Si II A6355) uniquely specify a SN la. (b) The spectrum and morphology of the 
host galaxy can identify it as an elliptical or SO, indicating that the supernovae is a type la, since 
only SNe la are found in these galaxy types. (Of course, the converse does not hold, since SNe 
la are also found in late spirals, at an even higher rate than in ellipticals, locally.) We will here 
consider E or SO host galaxies to indicate a SN la, with the caveat that it is possible, in principle, 
that someday a SN II may be found in these galaxy types, (c) The light curve shape can narrow 
the range of possible identifications by ruling out the plateau light curves of SNe IIP. (d) The 
statistics of the other classified supernovae discovered in the same search provide a probability 
that a random unclassified member of the sample is a SN la (given similar magnitudes above the 
detection threshold). 

Five of the seven supernovae discussed in this paper can be classified according to the criteria 
(a) and (b). Two are confirmed to be SNe la and three are consistent with SNe la: The spectrum 
of SN 1994an exhibits the major spectral features from 3700A to 6500A(SN rest frame wavelength) 
characteristic of SNe la ~3 days past maximum (rest frame time), including the Si II absorption 
near 6150A. SN 1994am was discovered in an elliptical galaxy, as identified by its morphology in 
a Hubble Space Telescope WFPC2 image (Figure ||), and by its spectrum, which matches that of 
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a typical present-day E or SO galaxy. We take these two SNe, 1994an and 1994am, to be SNe la. 
The spectrum of the SN 1994al host galaxy is a similarly good match to an E or SO, but until the 
morphology can be confirmed we consider SN 1994al to be consistent with a SN la. 

Spectra of both SN 1994F and SN 1994G are more consistent with an SN la spectrum at the 
appropriate number of rest-frame days (i S pect of Table |TJ) past the light curve maximum than any 
alternatives. The spectrum of SN 1994G was observed ~13 days past B maximum (rest frame) 
at the MMT by P. Challis, A. Riess, and R. Kirshner and ~15 days past maximum at the Keck 
telescope. The strengths of the Ca II H and K, Fe, and Mg II features closely match those of 
a SN la. The spectrum of SN 1994F was observed ~2 days before maximum (rest frame) by 
J. B. Oke, J. Cohen, and T. Bida during the commissioning of LRIS at the Keck telescope, and 
therefore was neither calibrated nor optimally sky-subtracted. The stronger SN I features (e.g. 
Ca II) do appear to be present nonetheless, so SN 1994F is unlikely to be a luminous SN II near 
maximum. The spectra for all of these supernovae will be reanalyzed and published once the 
late-time host-galaxy spectra are available, since host-galaxy features can confuse details of the 
supernova spectrum. 

To address the remaining two unclassified supernovae, we consider the classification statistics 
of our entire sample of >28 supernovae discovered to date by the Supernova Cosmology Project. 
Eighteen supernovae have been observed spectroscopically, primarily at the Keck 10 m telescope. 
Of these, the 16 supernovae at redshifts z > 0.35 discovered by our standard search technique are 
all consistent with a type la identification, while the two closer (z < 0.35) supernovae are type 
II. This suggests that >94% of the supernovae discovered by this search technique at redshifts 
z > 0.35 are type la. Moreover, the shapes of the light curves we observe are inconsistent with the 
"plateau" light curves of SNe IIP, further reducing the probability of a non-type la. In this paper, 
we will therefore assume an SN la identification for the remaining two unclassified supernovae. 

4.2. Photometry Reduction 

The two primary stages of the photometry reduction are the measurement of the supernova 
flux on each observed point of the light curve, and the fitting of these data to a SN la template 
to obtain peak magnitude and width (stretch factor). At high redshifts, a significant fraction 
(~50%) of the light from the supernova's host galaxy usually lies within the seeing-disk of the 
supernova, when observed from ground-based telescopes. We subtract this host galaxy light from 
each photometry point on the light curve. The amount of host galaxy light to be subtracted 
from a given night's observation is found from an aperture matched to the size and shape of 
the point-spread-function for that observation, and measured on the late-time images that are 
observed after the supernova has faded. 

The transmission ratio between the late-time image and each of the other images on the light 
curve is calculated from the objects neighboring the supernova's host galaxy that share a similar 
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Fig. 1. — Enlarged subimage showing the host galaxy of SN 1994am, and three neighboring 
galaxies before, during, and after the supernova explosion, (a) Isaac Newton 2.5-m Harris R- 
band observation on 19 December 1993, with 0.6 arcsec/pixel. (b) Kitt Peak 4-m Harris -R-band 
observation on 4 February 1994, with 0.47 arcsec/pixel, near maximum light of the supernova 
(indicated by arrow), (c) Cerro Tololo 4-m Harris i?-band observation on 15 October 1995, with 
0.47 arcsec/pixel. (d) Hubble Space Telescope (WFPC2) F814W /-band observation on 13 January 
1995, with 0.1 arcsec/pixel. At this resolution, the host galaxy of SN 1994am can be identified as 
an elliptical, providing strong evidence that SN 1994am is a type la supernova. 
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color (typically >25 objects are used). This provides a ratio that is suited to the subtraction of 
the host galaxy light. Another ratio between the images is calculated for the supernova itself, 
taking into account the difference between the color of the host galaxy and the color of the 
supernova at its particular time in the light curve by integrating host-galaxy spectra and template 
supernova spectra over the filter and detector response functions. (These color corrections were 
not always necessary since the filter-and-detector response function often matched for different 
observations.) The magnitudes are thus all referred to the late-time image, for which we observe 
a series of Landolt (1992) standard fields and globular cluster tertiary standard fields (Christian 
et al. 1985; L. Davis, private communication). The instrumental color corrections account for the 
small differences between the instruments used, and between the instruments and the Landolt 
standard filter curves. 

This procedure is checked for each image by similarly subtracting the light in apertures 
on numerous (~50) neighboring galaxies of angular size, brightness, and color similar to the 
host galaxy. If these test apertures, which contain no supernova, show flat light curves then the 
subtraction is perfect. The RMS deviations from flat test-lightcurves provide an approximate 
measurement, Otest; of photometry uncertainties due to the matching procedure together with 
the expected RMS scatter due to photon noise. We generally reject images for which these 
matching-plus-photon-noise errors are more than 20% above the estimated photon-noise error 
alone (see below). (Planned HST light curve measurements will make these steps and checks 
practically unnecessary, since the HST point-spread-function is quite stable over time, and small 
enough that the host galaxy will not contribute substantial light to the SN measurement.) 

Photometry Point Error Budget. We track the sources of photometric uncertainty at each step 
of this analysis to construct an error budget for each photometry point of each supernova. The 
test-lightcurve error, a t est, then provides a check of almost all of these uncertainties combined. The 
dominant source of photometry error in this budget is the Poisson fluctuation of sky background 
light, within the photometry apertures, from the mean-neighborhood-sky level. A typical mean 
sky level on these images is iV s k y = 10, 000 photoelectrons (p.e.) per pixel, and it is measured to 
approximately 2% precision from the neighboring region on the image. Within a 25-pixel aperture, 
the sky light contributes a Poisson noise of (25A r s k y ) 1 ^ 2 ~ 500 p.e. The light in the aperture from 
the supernova itself and its host galaxy is typically ~6,000 p.e., and contributes only ~77 p.e. 
of noise in quadrature, negligible compared to the sky background noise. Similarly, the noise 
contribution from the sky light in the subtracted late-time images after the supernova fades is 
also negligible, since the late time images are typically >9 times longer exposures then the other 
images. The supernova photometry points thus typically each have <7 p hoton ~ H% photon-noise 
measurement error. 

Much smaller contributions to the error budget are added by the previously-mentioned 
calibration and correction steps: The magnitude calibrations have uncertainties between 1% and 
4%, and uncertainties from instrumental color corrections are <1%. Uncertainties introduced in the 
"flat-fielding" correction for pixel-to-pixel variation in quantum efficiency are less than l/\/T0 of 
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the sky noise per pixel, since >10 images are used to calculate the "flat field" for quantum efficiency 
correction. This flat-fielding error leads to ~32 p.e. uncertainty in the supernova-and-host-galaxy 
light, which is an additional 5% uncertainty in quadrature above <7 p hoton- F° r most images the 
quadrature sum of all of these error contributions agrees with the overall test dispersion, <7 test . 
The exceptions mentioned above generally arise from point-spread-function variations over a given 
image in combination with extremely poor seeing, and such images are generally rejected unless 
otest i s within 20% of the expected quadrature sum of the error contributions. 

Since the target fields for the supernova search were chosen at high Galactic latitudes 
whenever possible, the uncertainties in our Galaxy extinction, An (based on values of E(B — V) 
from Burstein & Heiles 1982), are generally <1% on the photometry measurement. The one major 
exception is SN 1994al, for which there is more substantial Galactic extinction, Ar = 0.23 mag, 
and hence we quote a more conservative extinction error of ±0.11 mag for this supernova. 

All of these sources of uncertainties are included in the error budget of Table § The 
photometry data points derived from this first stage of the reduction are shown in Figure ^. 
Although we observe >2 images for each night's data point to correct for cosmic rays and pixel 
flaws, these images have been combined in producing these plots (and in the preceding noise 
calculations, for direct comparison). Since the error bars depend on the uncertainty in the 
host-galaxy measurement, there is significant correlation between them, particularly for nights 
with similar seeing. Therefore, the usual point-to-point statistics for uncorrelated data do not 
apply for these plots. Further details concerning the calibration observations, color corrections, 
and data reduction, with mention of specific nights, telescopes, and supernovae, are given in 
Perlmutter et al. (1995), and the forthcoming data catalog paper. 



4.3. Light Curve Fit 

At z ~ 0.4, the light that leaves the supernova in the B band arrives at our telescopes 
approximately in the R band. The second stage of the photometry reduction, fitting the light 
curve, must therefore be performed using a .fT-corrected SN la light curve template: We use 
spectra of several low-redshift supernovae, well-observed over the course of their light curves, to 
calculate a table of cross-filter K corrections, KBn(t), as a function of light curve time. These 
corrections account for the mismatch between the redshifted B band and the R band (including 
the stretch of the transmission- function width), as well as for the difference in the defined zero 
points of the two magnitude systems. For each high-redshift supernova, we can then construct a 
predicted i?-band template light curve as would be observed for the redshift, 7^ bserve ) based on 
the -B-band "standard" template (i.e., with Amis = 1-1 mag) as actually observed in the SN rest 
frame: 

^bserve (^ = Restframe ^ + Kjjr ^ (?) 
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Fig. 2. — ( Continued on next page.) i?-band light curve photometry points for the first seven high- 
redshift supernovae discovered by the Supernova Cosmology Project, and best-fit template SN la 
light curves. The left panels show the relative flux as function of observed time (not supernova- 
rest- frame time). The right panels show observed i?-band magnitude versus observed time. Note 
that there is significant correlation between the error bars shown, particularly for observations with 
similar seeing, since the error bars depend on the uncertainty in the host-galaxy measurement that 
have been subtracted from these measurements (see text). An /-band light curve is also shown for 
SN 1994G; other photometry points in / and B for the seven supernovae are not shown on this plot. 
The rising slope (in mag/day) of the template light curve before rest-frame day —10 (indicated by 
the grey part of the curves) is not well-determined, since few low-redshift supernovae are discovered 
this soon before maximum light. A range of possible rise times was therefore explored (see text). 




Fig. 2. — continued. 
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where the observed time dependence, t' = t(l + z), accounts for the time dilation of events at 
redshift z. The calculations of KsR(t) are described and tabulated, in Kim, Goobar, & Perlmutter 
(1996), with an error analysis that yields uncertainties of <0.04 mag for redshifts z < 0.6. Table 1 
lists the values of Kbr at maximum light for the redshifts of the seven supernovae. 

For the exceptional case of SN 1994al, with significant E{B—V) from our own Galaxy, we 
also calculate the extinction, Aft, as a function of supernova-rest-frame time, again using a series 
of redshifted supernova spectra multiplied by a reddening curve for an E(B— U)b&h value given 
by Burstein & Heiles (1982). For this particular supernova's redshift, An = 0.23 ± 0.11, with 
variations of 0.01 magnitude for the dates observed on the light curve. For the redshifts of all 
seven supernovae, we find essentially the same ratio TZr = Ar/ 'E(B—V)b&h = 2.58 ± 0.01, at 
maximum light. 

The photometry data points for each supernova are fit to the as-observed -R-band template, 
T^ bscrvc , with free parameters for the stretch-factor, s, the R magnitude at peak, uir, the date of 
peak, tpeaki and the additive constant flux, <? re sid) that accounts for the residual host-galaxy light 
due to Poisson error in the late-time image: 

/(*! m R , s, t peak , 5resid ) = l0-°-^-<+^^-W)] + , resid , ( 8 ) 

where m R is the calibration zero-point for the R observations and T^ bserve (t) is normalized to zero 
at t = 0. (For one of the supernovae with sufficient I band data, SN 1994G, the fitting function / 
also includes the redshifted, X-corrected U-band template, and an additional parameter for the 
R—I color at maximum is fit.) We fit to flux measurements, rather than magnitudes, because 
the error bars are symmetric in flux, and because the data points have the late-time galaxy light 
subtracted out and hence can be negative. In this fit, particular care is taken in constructing the 
covariance matrix (see, e.g., Barnett et al. 1996) to account for the correlated photometry error 
due to the fact that the same late-time images of the host galaxy are used for all points on a light 
curve. 

Light Curve Error Budget. We compute uncertainties for tur and s using both a Monte 
Carlo study and a mapping of the x 2 function; both methods yield similar results. Table [l] lists 
the best-fit values and uncertainties for thr and s. Usually, many data points contribute to 
the template fit, so the uncertainty in the peak, uir, is less than the typical ~11% photometry 
uncertainty in each individual point. However, our Monte Carlo studies show that it is usually 
very important that high quality late-time and pre-maximum data points be available to constrain 
both s and mpt- We have explored the consequences of adding or subtracting observations on the 
light curves, although the error bars reflect the actual light curve sampling observed for a given 
supernova. 

The rising slope of the template light curve before rest-frame day —10 (indicated by the 
grey part of the curves in Figure ^) is not well-determined, since few low-redshift supernovae 
are discovered this soon before maximum light. A range of possible rise times was therefore also 
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explored. Only two of the supernovae show any sensitivity to the choice of rise time. The effect is 
well within the error bars of the stretch factor, and negligible for the other parameters of the fit. 

For the analyses of this paper, we translate these observed R magnitudes back to the 
"effective" B magnitudes, mg = niR — Ar — Kbr, where all the quantities (see Table |l|) are for 
the light curve l?-band peak (SN rest frame), and we have corrected for Galactic extinction, Ar, 
at this stage. In the following section, we directly compare these "effective" B magnitudes with 
the B magnitudes of the low-redshift supernovae using the equations of Section 2, substituting 
the .ff-corrected effective B magnitudes, thb, and the B-h&nd zero point, Mb, for the bolometric 
magnitudes, m and Ai. 



5. Results for the High-Redshift Supernovae 

5.1. Dispersion and Width-Brightness Relation 

Before using a width-luminosity correction, it is important to test that it applies at high 
redshifts, and that the magnitude dispersions with and without this correction are consistent 
with those of the low-redshift supernovae. We study the peak-magnitude dispersion of the seven 
high-redshift supernovae by calculating their absolute magnitudes for an arbitrary choice of 
cosmology. This allows the relative magnitudes of the supernovae at somewhat different redshifts 
(z = 0.35-0.46) to be compared. The slight dependence on the choice of cosmology is negligible for 
this purpose, for a wide range of Qm and Choosing = 1 and £2^ = 0, the RMS dispersion 
about the mean absolute magnitude for the seven supernovae is om b = 0-27. (We find the same 
RMS dispersion for the best-fit cosmology discussed below in Section 5.3.) 

Figure || shows the difference between the measured and the "theoretical" (Equation ^| for 
= 1, = and M. = jMg^) uncorrected effective B magnitudes, mg — m^ eory , as a 
function of the best-fit stretch factor s or, equivalently, Amis. The zero-point intercept, M.g '^ rr , 
from Equation [| was used in the calculation of m^ eory , so that the data points could be compared 
to the Equation || width-luminosity relation found for the nearby sample of Hamuy et al. (the 
solid line of Figure |||). A different choice of (Om> ^a) or of the magnitude zero-point M. would 
move the line in the vertical direction relative to the points, so the comparison should be made 
with the shape of the relation, not the exact fit. The width-luminosity relation seen for nearby 
supernovae appears qualitatively to hold for the high-redshift supernovae. 

To make this result quantitative, we use the prescription of Equation [| to "correct" the 
magnitudes to that of a supernova with Amis = 1-1 mag, by adding magnitude corrections, 
^corV = — 0.86(Ami5 — 1.1), where we derive Amis from s using Equation ^. Note that the 
resulting corrected magnitude, thb, corr — f^B ~\~ ^cotti listed in Table [l], has an uncertainty that 
is not necessarily the quadratic sum of the uncertainties for m b and A^V > because the best-fit 
values for s and itir (and hence tub) are correlated in the light-curve template fit. 
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Fig. 3. — The difference between the measured and "theoretical" (Equation |I| for = 1, ^a = 0, 
and Ai = Mg 1 '^ ) B magnitudes (uncorrected for the width-luminosity relation) versus the best- 
fit stretch factor, s, for the high-redshift supernovae. The stretch factor is fit in the supernova 
rest frame, i.e., after correcting for the cosmological time dilation calculated from the host-galaxy 
redshift (see Equations |7| and ||). If different values of (Om, ^a) had been chosen the labels 
would change, but the data points would not vary significantly within their error bars, since the 
range of redshifts is not large for these supernovae. The upper axis gives the equivalent values of 
Amis = 1.96(s _1 — 1) + 1.07 (Eq. ^). The solid line shows the width-luminosity relation (Eq. |5|) 
found by Hamuy et al. (1995, 1996) for an independent set of 18 nearby {z < 0.1) SNe la, for which 
0.8 < Amis < 1.75 mag. This range of light-curve widths is indicated by the shaded region. The 
curve and data points outside of this range are plotted in a different shade to emphasize that the 
relation is only established within this range. A different choice of (Om, or of the magnitude 
zero-point M would move the line in the vertical direction relative to the points. 
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Two of the high-redshift supernovae (shown with different-shaded symbols in Figure || 
and subsequent figures) have values of Amis that are outside of the range of values studied in 
the low-redshift supernova data set of Hamuy et al. The width-luminosity correction for these 
supernovae is therefore less reliable than for the other five supernovae, since it depends on an 
extrapolation. It is also possible that the true values of Amis for the two supernovae fall within 
this range, given the fit error bars. (One of these, SN 1992bi at z = 0.458, has very large error 
bars because its light curve has larger photometric uncertainties and the lightcurve sampling was 
not optimal to constrain the stretch factor s.) In Table [l], the corrections, A^$, and corrected 
magnitudes, mB )C om for these supernovae are therefore listed in brackets, and we also give the 
corrected magnitudes, m^* 1 ^™ 6 , obtained using only the most extreme A^^ corrections found for 
low-redshift supernovae. In the following analysis we calculate all results using just the other five 
supernovae that are within the Hamuy et al. range of Amis. As a cross-check, we then provide 
the result using all seven supernovae, but without width-luminosity correction. 

The peak-magnitude RMS dispersion for the five high-redshift supernovae that are within 
the Amis range improves from <jm b = 0.25 mag before applying the width-luminosity correction 
to (Tm b ,cott = 0.19 mag after applying the correction. These values agree with those for the 18 
low-redshift Calan/Tololo supernovae, c^ nuy = 0.26 and (T^^orr = 0-17 m ag> suggesting that the 
supernovae at high-redshift are drawn from a similar population, with a similar width-luminosity 
correlation. Note that although the slope of the width-luminosity relation, Equation [|, stays the 
same, it is possible that the intercept could change at high redshift, but this would be a remarkable 
conspiracy of physical effects, since the time scale of the event appears naturally correlated with 
the strength and temperature of the explosion (see Nugent et al. 1996). 



5.2. Color and Spectral Indicators of Intrinsic Brightness 

With the supplementary color information in the / and B bands, and the spectral information 
for three of the supernovae, we can also begin to test the other indicators of intrinsic brightness 
within the SN la family. The broadest, most luminous supernova of our five-supernova subsample, 
SN 1994H, is bluer than B—R = 0.7 at the 95% confidence level 2 ± 2 days before maximum light 
(this is only a limit because of possible clouds on the night of the B calibration). For comparison, 
Nugent (1996, private communication) has estimated observed B—R colors within ~4 days of 
maximum light as a function of redshift using spectra for a range of SN la sub-types, from the 
broad, superluminous SN 1991T to the narrow, subluminous SN 1991bg. SN 1994H is only slightly 
bluer than the B—R ~ 1.0 ± 0.15 color of SN 1991T at z = 0.374 (where the error represents the 
range of possible host-galaxy reddening). However, SN 1994H does not agree with the redder 
colors (at z = 0.374) of SN 1981B (B-R « 1.7 mag) and SN 1991bg (B-R ra 2.6 mag). SN1991T 
has a broad light curve width of Amis = 0.94 ± 0.07, which agrees within error with the light 
curve width of SN 1994H, Ami 5 = 0.91 ± 0.09, whereas SN 1981B and SN 1991bg both have 
"normal" or narrower light curve widths, Amis = 1-10 ± 0.07 and Amis = 1-93 ± 0.10 mag. 
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The multi-band version of the "correction template" analysis (Riess et al. 1996) is designed 
to fit simultaneously for the host galaxy extinction (discussed later) and the intrinsic brightness of 
the supernova. The rest frame V light curve (approximately redshifted into our observed I) is the 
strongest indicator of the SN la family parameterization when using this approach. We thus use 
this technique to analyze SN1994G, with its well sampled / light curve in addition to R data. We 
find that the supernova is best fit by a correction template that indicates that it is intrinsically 
overluminous by 0.06 ± 0.38 magnitudes, compared to a Leibundgut-template supernova. This is 
consistent with the 0.05 ± 0.22 overluminosity found using the width-luminosity correlation and 
the best-fit stretch factor. 

(The "correction template" fit gives larger error bars than a simple stretch-factor fit because 
the correction templates have significant uncertainties, with day-to-day correlations that are, 
unfortunately, not tabulated. Currently, the fits to £?-band correction templates have even 
larger uncertainties than T/-band, possibly due to a poor fit of the linear correction model to the 
low-redshift data. This approach is therefore not useful for the other six high-redshift supernovae 
of this first set, since these large correction-template uncertainties propagate into >1 magnitude 
uncertainties in the fits to the observed R band data.) 

Both intrinsically fainter supernovae and supernovae with host galaxy extinction can appear 
redder in R—I (corresponding to approximately B—V in the supernova rest frame). For the 
several supernova for which we have scattered /-band photometry, we thus can use it to estimate 
extinction only after the supernova sub-type has been determined using the stretch factor. 
SN 1994G fit closely to the s = 1 template. An s = 1 supernova at z = 0.420 will have an expected 
observed color R — I = 0.16 ±0.05 mag at observed- R maximum light. For SN 1994G, we observed 
R-I = -0.12 ± 0.17 mag. This gives E(R-I) ohseived « E(B - F) rcstframe = -0.28 ± 0.18 mag, 
or Ay < 0.01 mag at the 95% confidence level. The Riess et al. "correction template" analysis of 
SN 1994G also yields a bound on extinction of Ay < 0.01 mag (90% confidence). 

The more recent high-redshift supernovae studied by the Supernovae Cosmology Project have 
more complete color and spectral data, and so should soon provide further tests of these additional 
SN la luminosity indicators. In addition, the spectrum of SN 1994an covers both wavelength 
regions in which the Nugent et al. (1996) line ratios correlate well with supernova magnitudes. 
This analysis awaits the availability of the spectrum of the host galaxy without the supernova 
light, so that we can subtract the galaxy-spectrum "contamination" from the supernova spectrum. 
Further late-time B-b&nd images of SN 1994an will also allow B—R color cross-checks using B 
images that were observed four days after maximum. 

5.3. Magnitude-Redshift Relation and the Cosmological Parameters 

Figure 0(a) shows the Hubble diagram, tub versus logcz, for the seven high-redshift 
supernovae, along with low-redshift supernovae of Hamuy et al. (1995) for visual comparison. The 
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solid curves are plots of m^ eory , i.e., Equation [| with M. = M.b for three A = cosmologies, 
(f^Mj ^a) = (0) 0), (1, 0), and (2, 0); the dotted curves, which are practically indistinguishable 
from the solid curves, are for three flat cosmologies, (Ami ^a) = (0.5, 0.5), (1, 0), and (1.5, —0.5). 
These curves show that the redshift range of the present supernova sample begins to distinguish 
the values of the cosmological parameters. Figure |](b) shows the same magnitude-redshift relation 
for the data after adding the width-luminosity correction term, A^'; the curves of m^ eory in (b) 
are calculated for M. = ■Mb^tt- 

Figure | shows the 68% (la), 90%, and 95% (2a) confidence regions on the Qm _ ^A plane for 
the fit of Equation |J to the high-redshift supernova magnitudes, mB iCorn after width-luminosity 
correction, using M. = M.% 1 '^ rr . In this fit, the magnitude dispersion, cr^^^ rr = 0.17, is added 
in quadrature to the error bars of mB )Cori to account for the residual intrinsic dispersion, after 
width-correction, of our model (low-redshift) standard candles. 

The two special cases represented by the solid lines of Figure || yield significant measurements: 
We again fit Equation |I| to the high-redshift supernova magnitudes, fnB iCorr , but this time with 
the fit constrained first to a A = cosmology (the horizontal line of Figure ||), and then to a 
flat universe (the diagonal line of Figure ||, with Sltotai = + ^A = 1)- In the case of a A = 
cosmology, we find the mass density of the universe to be S7m = 0.88 IqJo- ^ or a unrverse 
(the diagonal line of Figure ||, with fitotai = + ^A = 1), we find the cosmological constant to 
be f^A = 0.06 ^0 34, consistent with no cosmological constant. For comparison with other results 
in the literature, we can also express this fit as a 95% confidence level upper limit of < 0.51. 
Finally, this fit can be described by the value of the mass density, J7m = 1 — ^A = 0.94 <28- 
(For brevity, we will henceforth only quote the Q\ fit results for the flat-universe case, since either 
Q\ or can be used to parameterize the flat-universe fit.) The goodness-of-fit is quantified by 
Q(x 2 W), the probability of obtaining the best-fit x 2 ° r higher for v = A^Ne — 1 degrees of freedom 
(following the notation of Press et al. 1986, p. 165). For both the A = and flat universe cases, 
Q( X 2 W) = 0.75. 

The error bars on the measurement of f^A or Qm for a flat universe are about two times 
smaller than the error bars on S7m for a A = universe. This can be seen graphically in Figure ||, 
in which the confidence region strip makes a shallow angle with respect to the f^A = line but 
crosses the f2 tota i = 1 line at a sharper angle. Note that these error bars in the constrained 
one-dimensional fits for A = or Ototai = 1 are smaller than the intersection of the constraint 
lines of Figure [| and the 68% contour band. This is because a different range of A% 2 = X 2 ~ Xmin 
corresponds to 68% confidence for one free parameter, v = 1 degree of freedom, as opposed to two 
free parameters, v = 2 (see Press et al. 1986, pp. 532-7). 

We also analyzed the data for the same five supernovae "as measured," that is, 
without correcting for the width-luminosity relation, and taking the uncorrected zero point, 
M B = -3.17 ± 0.03 of Hamuy et al. (1996). We find essentially the same results: Q M = 0.93 ±oH 
for A = and Q\ = 0.03 34 for a flat Ototai = 1 universe. The measurement uncertainty for 
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Fig. 4. — Hubble diagrams for the first seven high-redshift supernovae, (a) uncorrected rag, 
with low-redshift supernovae of Hamuy et al. (1995) for visual comparison; (b) mg ]H)rr after 
"correction" for width-luminosity relation. The square points are not used in the analysis, because 
they are corrected based on an extrapolation outside the range of light curve widths of low-redshift 
supernovae (see text and Table [l]). Insets show the high-redshift supernovae on magnified scales. 
The solid curves in (a) and (b) are theoretical m# for (Om, ^a) = (0, 0) on top, (1, 0) in middle, and 
(2, 0) on bottom. The dotted curves, which are practically indistinguishable from the solid curves, 
are for the flat universe case, with (Qmj ^a) = (0.5, 0.5) on top, (1, 0) in middle, and (1.5, —0.5) on 
bottom. The inner error bars on the data points show the photometry measurement uncertainty, 
while the outer error bars add the intrinsic dispersions found for low-redshift supernovae, cr^ muy 
for (a) and ^/^corr f° r O 3 )) f° r comparison to the theoretical curves. Note that the zero-point 
magnitude used for (b), A^g'^ rr , and hence the effective ms scale, is shifted slightly from the 
uncorrected Mb used for the curves of (a). 
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Fig. 5. — Contour plot of the 68% (la), 90%, and 95% (2d) confidence regions in the f^A versus 
r^M plane, for the first seven high-redshift supernovae. The solid lines show the two special cases 
of a A = cosmology and a flat (Sltotal = + = 1) cosmology. (Note that the constrained 
one-dimensional confidence intervals for A = or ^totai = 1 are smaller than the intersection of 
these lines with the two-dimensional contours, as discussed in the text.) The two labeled corners 
of the plot are ruled out because they imply: (upper left corner) a "bouncing" universe with no 
big bang (see Carroll et al. f992), or (lower right corner) a universe younger than the oldest heavy 
elements, to < 9.6 Gyr (Schramm 1990), for any value of Hq > 50 km s _1 Mpc -1 . 
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the analysis using light-curve-width correction is not significantly better than this one, without 
correction, because for this particular data set the smaller calibrated dispersion, (tm b ,cott as 
opposed to ctm b , is offset by the uncertainties in the measurements of the light curve widths. The 
goodness of fit, however, for the uncorrected magnitudes is lower, Qix^W) = 0.34, for both the 
A = and flat universe. This indicates that the width-luminosity relation provides a better model, 
although this was already seen in the improved RMS dispersion of the corrected magnitudes. 
Table § summarizes the results. 

As a cross-check, we also analyzed the results for all seven supernovae, not "width-corrected," 
including the two that had measured widths outside the range of the low-redshift width correction. 
We find n M = 0.70 for A = and ft A = 0.15 ±gj| for O tota i = 1. 

We emphasized in Section 2 that the la confidence region of Figure || is not parallel to the 
contours of constant qo. For comparison, the qo contours are drawn (dashed lines) in Figure ||. 
The closest qo contour varies from ~0.5 to ~1 in the region with Qm < 1- Any single value of qo 
would thus be only a rough approximation to the true confidence interval dependent on J7m and 
Oa shown in Figure |f| 

6. Checks for Sources of Systematic Error 

Since this is a first measurement of cosmological parameters using SNe la, we list here some 
of the more important concerns, together with the checks and tests that would address them. By 
considering each of the separate subsamples of SNe that avoid a particular source of bias, we can 
show that no one source of systematic error alone can be responsible for the measured values 
of and Q\. We cannot, however, exclude a conspiracy of biases each moving an individual 
supernova's measured values into agreement with this result. Fortunately, this sample of seven 
supernovae is only the beginning of a much larger data set of SNe la at high redshifts, with more 
complete multiband light curves and spectral coverage. 

Width- Brightness Correction. For the present data set, the measurements of f^M and Oa 
are the same whether or not we apply the light-curve width-luminosity correction, and so do not 
depend upon our confidence in this empirical calibration. This is due to the similar distribution 
of light curve widths for the high-redshift supernovae and the low-redshift supernovae used for 
calibration. 

For future data sets, it is possible that the width distribution will differ, for example if we were 
to find more supernovae in clusters of ellipticals and confirm the tendency to find narrower /fainter 
supernovae in ellipticals. (Note that SN 1994am, in an elliptical host, does in fact have a somewhat 
narrow light curve, with s = 0.84.) For such a data set the result would depend on the validity of 
the width-luminosity correction. Such a correction dependence could be checked by restricting the 
analysis to the supernovae that pass the Vaughan et al. (1995) B—V color test for "normals," and 
then applying no width-luminosity correction. 
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Fig. 6. — Contour plot of the la (68%) confidence region in the versus Qm plane, for the first 
seven high-redshift supernovae. The solid lines show contours for the age of the universe (in Gyr), 
normalized to Hq = 70 in units of km s _1 Mpc -1 . The dashed lines are the contours of constant 
deceleration parameter qo = 0.0 (top), 0.5 (middle), 1.0 (bottom). 
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Extinction. While the extinction due to our own Galaxy has been incorporated in this 
analysis, we do not know the SN host-galaxy extinctions. Note that correcting for any neglected 
extinction for the high-redshift supernovae would tend to brighten our estimated SN effective 
magnitudes and hence move the best fit of Equation |] towards even higher Qm an d lower Q\ 
than the current results. We can check that the extinction is not strongly affecting the results by 
considering two supernovae for which there is evidence against significant extinction. SN 1994am 
is in an elliptical galaxy, and for SN 1994G there is the previously-mentioned R—I color that 
provides evidence against significant reddening. The best fit values for this two-SN subset are 
consistent with that of the full set of SNe: Q M = 1.31 tlH for A = and A = -0.15 for 

^total = 1- 

If there were uncorrected host-galaxy extinction for the low-redshift supernovae used to 
find the magnitude zero point M, this would lead to an opposite bias. The 18 Calan/Tololo 
supernovae, however, all have unreddened colors (B — V < 0.2 mag). For the range of widths 
of these supernovae (0.8 < Amis ^ 1-75 mag), the range of intrinsic color at maximum light 
should be —0.05 < B — V < 0.5 mag, so the color excess is limited to E{B — V) < 0.25 mag. A 
stronger constraint can be stated for the seven of these 18 Calan/Tololo supernovae that were 
included in a sample studied by Riess et al. (1996), who used the multicolor correction-template 
method to estimate extinction. They found that only one of these seven supernovae showed any 
significant extinction. (For that one supernova, SN 1992P, they found Ay = 0.11 mag.) It is thus 
unlikely that host-galaxy extinction of the Calan/Tololo supernovae is strongly biasing our results, 
although it will be useful to test the remaining supernovae of this set for evidence of extinction, 
using the multicolor correction-template method. (It should be noted that the method of Riess 
et al. estimates total extinction values for several supernovae that are significantly less than the 
Galactic extinction estimated by Burstein & Heiles 1982, hence some caution is still necessary in 
interpreting these results.) 

A general argument can be made even without these color measurements: the intrinsic 
magnitude dispersions, ^^corr = 0-17 or our value o"A/ s ,corr = 0.19 mag, provide an upper bound 
on the typical extinction present in either the low redshift or the high redshift supernova samples, 
since a broad range of host-galaxy extinction (which would have a larger mean extinction) would 
inflate these dispersions. We use Monte Carlo studies to bound the amount of bias due to a 
distribution of extinction that would inflate a hypothetical intrinsic dispersion of a ^yp° tliesis = q.13 
to the values actually measured, o^^cott = 0-17 ^ anc ^ a M B ,corr = 0.19 ± 0.08 mag. We find, 
at the 90% confidence level, less than 0.06 magnitude bias towards higher Qm arid lower Q\ and 
less than 0.10 magnitude bias towards lower S7m and higher These measurements and bounds 
on extinction will be even more important as we study supernovae at still higher redshifts and 
need to check for evolution in the host galaxy dust. 

We choose high Galactic latitude fields whenever possible, so that the extinction correction 
for our own Galaxy and its uncertainty are negligible. The one major exception in this current 
data set is SN 1994al, with Ar ~ 0.23 mag. (This supernova also appears to be a somewhat 
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fainter outlier among the width-corrected magnitudes on Figure |4]b.) Since there is more 
uncertainty in the Galactic extinction correction for this supernova, we have fit just the other four 
width-corrected supernova magnitudes. For SN 1994H, SN 1994am, SN 1994G, and SN 1994an, 
we find U M = 1.25 ±{j;j$ for A = and n A = -0.12 t H for O total = 1. 

K Correction. The generalized K correction used to transform i?-band magnitudes of high 
redshift SNe to l?-band rest frame magnitudes was tested for a variety of SN la spectra and found 
to vary by less than 0.04 mag for redshifts z < 0.6 (Kim, Goobar, & Perlmutter 1996). The 
test spectra sample included examples in the range of light-curve widths between Amis = 0.94 
mag (for SN 1991T) and Amis = 1-47 mag (for SN 1992A), but K corrections still need to be 
calculated to check for possible difference outside of this regime of the SN la family. Two of the 
five supernovae we used in our measurement fall just outside of this range of widths. As a check 
of a possible systematic error due to this, we consider just the subsample of three supernovae, 
SN 1994al, SN 1994am, and SN 1994G, with light curve widths within the studied range of K 
corrections. We find M = 0.61 toil for A = and ft A = 0.18 l£ff for O total = 1. 

Malmquist Bias. The tendency to find the most luminous members of a distribution at 
large distances in a magnitude-limited search would appear to bias our results towards larger 
values of Om and lower values of Qa- However, as the SN la Hubble diagram uses the peak 
brightness rather than that at detection, and the supernovae are "corrected" from different 
intrinsic brightnesses, any such Malmquist bias would not operate the same way it does for a 
population of normally-distributed static standard candles. The key issue is whether supernova 
samples are strongly biased at the detection level, as would be the case if most supernovae were 
discovered close to the threshold value. Our detection-efficiency studies (see Pain et al. 1996) show 
that the seven high-redshift supernovae were detected approximately 0.5 to 2 magnitudes brighter 
than our limiting (50% efficiency) detection magnitude, m d _ 50 (see the final row of Table [l]), and 
the efficiency on our CCD images drops off slowly beyond this limit. (Several of the SNe that 
we find at a significantly brighter magnitude than our threshold are in clusters that are closer 
than the limiting distance for a SN la. Such inhomogeneities in galaxy distribution can lead to a 
sample of supernovae that are not distributed primarily near the magnitude limit as expected in a 
magnitude limited search.) In contrast, the Calan/Tololo survey detected most of the low-redshift 
supernovae within ~0.7 magnitude of detection threshold and their efficiency on photographic 
plates dropped quickly beyond that limit (Maza, private communication). Malmquist bias may 
therefore result, counter- intuitively, in a slightly more luminous sample of the intrinsic SN la 
distribution for the low-redshift photographic search than for the high-redshift CCD search. 

Since, however, our current results already suggest a relatively high value for Qm and 
low value for Oa, we have nonetheless checked the possibility of Malmquist bias distorting the 
distributions of magnitudes we find at high redshift. First, as a rough test, we have compared 
the results (uncorrected for light-curve width) for the three supernovae discovered closest to the 
50%-efficiency detection threshold, m d 50 , to the three supernovae discovered farthest from the 
threshold. We find that the difference in measured J7m and Oa between these two subsets is not 
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significant, and it is opposite to the direction of Malmquist bias. 

A more detailed quantitative study was based on a Monte Carlo analysis, in which the 
detection efficiency curve for each supernova was used with the Calan/Tololo "corrected" 
magnitude dispersion, (^m^ovv ^° estimate the magnitude bias that should be present for each of 
the seven fields on which we discovered high-redshift supernovae. We find that only one supernova, 
SN 1994G, is on a field that shows any significant magnitude bias. Even for this supernova the 
bias, 0.01 mag, is still well below the intrinsic dispersion of the supernovae. Hence, when we 
reanalyze the data set, correcting SN 1994G by this amount, we find an insignificant change in our 
results for J7m and 

Because the high-redshift supernovae include both intrinsically narrow, subluminous cases 
and intrinsically broad, overluminous cases at comparable redshifts, a simple cross-check for 
Malmquist bias that is independent of detection-efficiency studies can be made by comparing the 
results for these two subsets separately. Malmquist bias would affect these two subsets differently, 
leading to a higher Qyi and lower Q,^ for the broad, overluminous subsample than for the narrow, 
subluminous subsample. With our current sample we can compare only two supernovae in each 
of these subsamples that are in the "correctable" range of Amis, so this will be a particularly 
interesting test to apply to the full sample of high-redshift supernovae when their light-curve 
observations are completed. The current data show no evidence of Malmquist bias. 

Supernova Evolution. Although there are theoretical reasons to believe that the physics of 
the supernova explosion should not depend strongly on the evolution of the progenitor and its 
environment, the empirical data are the final arbiters. Both the low-redshift and high-redshift 
supernovae are discovered in a variety of host galaxy types, with a range of histories. The 
small dispersion in intrinsic magnitude across this range, particularly after the width-luminosity 
correction, is itself an indication that any evolution is not changing the relationship between the 
light curve width/shape and its absolute brightness. (Note that the one supernova in an identified 
elliptical galaxy gives results for the cosmological parameters consistent with the full sample of 
supernovae; such a comparative test will of course be more useful with the larger samples.) In 
the near future, we will be able to look directly for signs of evolution in the >18 spectra observed 
for the larger sample of high-redshift supernovae. So far, the spectral features studied match the 
low-redshift supernova spectra for the appropriate day on the light curve (in the supernova rest 
frame), showing no evidence for evolution. A more detailed analysis will soon be possible, as the 
host galaxy spectra are observed after the supernovae fade, and it becomes possible to study the 
supernova spectra without galaxy contamination. 

Gravitational Lensing. Since the mass of the universe is not homogeneously distributed, 
there is a potential source of increased magnitude dispersion, or even a magnitude shift, due to 
overdensities (or underdensities) acting as gravitational lenses that amplify (or deamplify) the 
supernova light. This effect was analyzed in a simplified "swiss cheese" model by Kantowski, 
Vaughan, & Branch (1995), and more recently using a perturbed Friedmann-Lemaitre cosmology 
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by Frieman (1996) and an n-body simulation of a ACDM cosmology by Wambsganss et al. (1996). 
The conclusion of the most recent analyses is that the additional dispersion is negligible at the 
redshifts z < 0.5 considered in this paper: Frieman estimates an upper limit of less than 0.04 
magnitudes in additional dispersion. Any systematic shift in magnitude distribution is similarly 
small: Wambsganss et al. give the difference between the median of the distribution and the true 
value for their particular mass-density distribution, which corresponds to a ~0.025 magnitude 
shift at z = 0.5. We can take this as a bound on the magnitude shift, since our averaged results 
should be closer to the true value than the median would be. Alternative models for the mass 
density distribution must satisfy the same observational constraints on dark matter power at small 
scales (from pairwise peculiar velocities and abundances of galaxy clusters) and therefore should 
give similar results. 

7. Discussion 

We wish to stress two important aspects of this measurement. First, although we have 
considered many potential sources of error and possible approaches to the analysis, the essential 
results that we find are independent of almost all of these complications: this direct measurement 
of the cosmological parameters can be derived from just the peak magnitudes of the supernovae, 
their redshifts, and the corresponding K corrections. 

Second, we emphasize that the high-redshift supernova data sets provide enough detailed 
information for each individual supernova that we can perform tests for many possible sources of 
systematic error by comparing results for supernova subsets that are affected differently by the 
potential source of error. This is a major benefit of these distance indicators. 

There have been many previous measurements and limits for the mean mass density of the 
Universe. The methods that sample the largest scales via peculiar velocities of galaxies and their 
production through potential fluctuations have tended to yield values of J7m close to unity (Dekel, 
Burstein, & White 1996 but see Davis et al. 1996 and Willick et al. 1997). Such techniques, 
however, only constrain f3 = Q^f jb where b is the biasing factor for the galaxy tracers, which is 
generally unknown (but see Zaroubi et al. 1996). Alternative methods based on cluster dynamics 
(Carlberg et al. 1995) and gravitational lensing (Squires et al. 1996) give consistently lower values 
of —0.2-0.4. There are many reasons, however, why these values could be underestimates, 
particularly if clusters are surrounded by extensive dark halos. The measurement using supernovae 
has the advantage of being a global measurement that includes all forms of matter, dark or 
visible, baryonic or non-baryonic, "clumped" or not. Our value of J7m = 0.88 ^o'go f° r a A = 
cosmology does not yet rule out many currently favored theories, but it does provide a first global 
measurement that is consistent with a standard Friedmann cosmology. Our measurement for a 
flat universe, = 0.94 Iq^I' ^ s more constraining: even at the 95% confidence limit, £Im > 0.49, 
this result requires non-baryonic dark matter. 
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In a flat universe (O tota i = 1) 5 the constraint on the cosmological constant from gravitational 
lens statistics is A < 0.66 at the 95% level (Kochanek 1995; see also Rix 1996 and Fukugita et al. 
1992). (In establishing this limit, Kochanek included the statistical uncertainties in the observed 
population of lenses, galaxies, quasars, and the uncertainties in the parameters relating galaxy 
luminosity to dynamical variables; he also presented an analysis of the various models of the 
lens, lens galaxy, and extinction to show the result to be insensitive to the sources of systematic 
error considered.) Our upper limit for the flat universe case, Oa < 0.51 ( 95% confidence level), 
is somewhat smaller than the gravitational lens upper limit, and provides an independent, direct 
route to this important result. The supernova magnitude-redshift approach actually provides a 
measurement, Oa = 0.06 1q 1s4' ra ther than a limit. The probability of gravitational lensing varies 
steeply with Oa down to the Oa ~ 0.6 level, but is less sensitive below this level, so the limits 
from gravitational lens statistics are unlikely to improve dramatically. The uncertainty on the Oa 
measurement from the high-redshift supernovae, however, is likely to narrow by \/3 as we reduce 
the data from the next ~18 high-redshift supernovae that we are now observing. 

For the more general case of a Friedmann-Lemaitre cosmology with the sum of Om and Oa 
unconstrained, the measured bounds on the cosmological constant are, of course, broader. A 
review of the observational constraints on the cosmological constant by Carroll, Press, & Turner 
(1992) concluded that the best bounds are —7 < Oa < 2, based on the existence of high-redshift 
objects, the ages of globular clusters and cosmic nuclear chronometry, galaxy counts as a function 
of redshift or apparent magnitude, dynamical tests (clustering and structure formation), quasar 
absorption line statistics, gravitational lens statistics, and the astrophysics of distant objects. 
Even in the case of unconstrained Om and Oa, we find a significantly tighter bound on the lower 
limit: Oa > —2.3 (at the two-tailed 95% confidence level, for comparison). Our upper limit is 
similar to the previous bounds, but provides an independent check: Oa < 11 for Om < 1 or 
A < 2.1 for M <2. 

These current measurements of Om and Oa are inconsistent with high values of the Hubble 
constant (Hq > 70 km s _1 Mpc -1 ) taken together with ages greater than to = 13 billion years for 
the globular cluster stars. This can be seen in Figure |6| which compares the confidence region for 
the high-redshift supernovae on the Om~Oa plane with the to contours. Above the 68% confidence 
region, only a small fraction of probability lies at ages 13 Gyr or higher for Hq > 70 km s _1 
Mpc -1 , and this age range favors low values for Om- Once again, the analysis of our next set of 
high-redshift supernovae will test and refine these results. 
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Table 1. Supernova Data and Photometry Error Budget 





SN 1992bi 


SN 1994H 


SN 1994al 


SN 1994F 


SN 1994am 


SN 1994G 


SN 1994an 


z [a] 


0.458 


0.374 


0.420 


0.354 


0.372 


0.425 


0.378 


m R 


22.01 (9) 6 


21.38 (5) 


22.42 (6) 


22.06 (17) 


21.73 (6) 


21.65 (16) 


22.02 (7) 


A R W 


0.003 (1) 


0.039 (4) 


0.228 (114) 


0.010 (1) 


0.039 (4) 


0.000 (1) 


0.132 (13) 


K BR 


-0.70 (1) 


-0.58 (3) 


-0.65 (1) 


-0.56 (3) 


-0.58 (2.5) 


-0.66 (1) 


-0.59 (2.5) 


m B 


22.71 (9) 


21.92 (6) 


22.84 (13) 


22.61 (18) 


22.27 (7) 


22.31 (16) 


22.48 (7) 


s 


1.45 (18) 


1.09 (5) 


0.95 (12) 


0.67 (15) 


0.86 (3) 


1.01 (13) 


0.77 (9) 


Amis 


0.47 (17) 


0.91 (9) 


1.17 (27) 


2.04 (65) 


1.39 (10) 


1.04 (25) 


1.65 (32) 


A {ii} 


[0.55 (20)} d 


0.16 (9) 


-0.06 (23) 


[-0.81 (59)] d 


-0.25 (10) 


0.05 (22) 


-0.47 (30) 


?71_B,corr 


[23.26 (24)] d 


22.08 (11) 


22.79 (27) 


[21.80 (69)] d 


22.02 (14) 


22.36 (35) 


22.01 (33) 


t spect w 


84 (3) 


219 (2) 


209 (2) / 


-2 (2) 9 


203 (2/ 


13 (l) h 


3(2) 


Bands 


7? 


B,R 


R 


R 


R 


B,R,I 


B, R 




1.0 


1.9 


0.8 


0.4 


0.9 


0.5 


1.5 



Note. — The uncertainties in the least significant digit are given in parentheses. Section 4 of the text defines the 
variables. 

a The redshifts were measured from host-galaxy spectral features, and their uncertainties are all <0.001. 

b This value for tur of SN 1992bi incorporates more recent light curve and calibration data than were available 
when Perlmutter et al. (1995) was prepared. These new data have improved the photometry error bars, as that 
paper suggested. The small change in the tur value quoted is also partly due to the extra degree of freedom in the 
new fit of the template light curve to a stretch factor, s. 

c Extinction for our Galaxy, Ar = 1ZrE(B—V)b&h, where TZr = 2.58 ± 0.01 was calculated using well-observed 
SN la spectra redshifted appropriately for each supernova, and the E(B— V)b&h values for each supernova coordinate 
are from Burstein & Heiles (1982) and Burstein, private communication (1996). 

d Since SN 1992bi and SN 1994F have Amis values outside of the range for low-redshift SNe, we do not use these 
lightcurve-width-corrected values for m_B iCor r for any of the results of this paper. In the figures, we do plot them as 
grey squares, along with white squares for the following alternative values that are corrected only to the extreme 
of A'orJ 1 values for low-redshift supernovae: = 22.97(24) for SN 1992bi, and m c fl xt c r ™ c = 22.05(69) for 

SN 1994F. 

e ts P ect is the number of days past B maximum in the supernova rest frame of the first spectrum. Generally, only 
the host galaxy spectrum is bright enough to be useful later than t spcc t > 30 days. 

f Spectrum of host galaxy showed the strong 4000A break of an elliptical galaxy, indicating that SN 1994al and 
SN 1994am are type la supernovae (see also Figure Q). 

g Keck spectrum of SN 1994F observed by J.B. Oke, J. Cohen, and T. Bida. 

h MMT spectrum observed by P. Challis, A. Riess, and R. Kirshner. We also observed a spectrum at the Keck 
telescope 15 rest-frame days past B maximum. 

1 m d _ 50 is the magnitude difference between the discovery magnitude and the detection threshold, i.e. the magnitude 
at which 50% of simulated supernovae are detected for the image in which the supernova was discovered. Note that 
the discovery magnitude is generally fainter than the peak magnitude, rriR, since most discoveries are before maximum 
light. 
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Table 2. Cosmological Parameters from SNe 94H, 94al, 94am, 94G, and 94an 





for no-A universe 
(Oa = 0) 


for fiat universe 

(Ototal = + = 1) 


Light-Curve- Width Corrected" 
Uncorrected 6 


n aa + - 69 

U - 88 -0.60 
n OQ +0 69 
u - ao -0.60 


0.06 ±° f 4 
OS +0 ' 29 

U.UO _o.34 



a The results are the best fit of the lightcurve- width-corrected data, ms iCorr = 
m B + A^rr 1 , of Table § to Equation §, with M = M^'cit = ~ 3 - 32 ± °- 05 - The 
probability of obtaining the best-fit x 2 or higher for both the A = and flat 
universes is Q(x 2 |^) — 0.75. 

b The results are the best fit of the data, ms, of Table [l] to Equation ^, with 
M — Mb = —3.17 ± 0.03. The probability of obtaining the best-fit x 2 or higher 
for both the A = and flat universes is Q(x 2 \ l/ ) = 0.34. 



